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Abstract Non-Markovian effects upon the Brownian movement of a free particle in 
the presence as well as in the absence of inertial force are investigated within the 
framework of Fokker — Planck equations (Rayleigh and Smoluchowski equations). 
More specifically, it is predicted that non-Markovian features can enhance the 
values of the mean square displacement and momentum, thereby assuring the 
mathematical property of differentiability of the these physically observable 
quantities. 
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1 Introduction 



The Brownian motion of a particle with mass m and position X = X(i) may be 
described by the Langevin equation [1] 

d 2 X dV(X) dX 

where the inertial force md 2 X/dt 2 offsets a set of three sorts of forces: a 
conservative force derived from an external potential V(X); a linearly velocity- 
dependent dissipative force, —ydX/dt, accounts for stopping the particle's motion; 
and a stochastic force L(t), dubbed Langevin's force, responsible for activating the 
particle's movement. The damping constant y exhibits dimensions of mass per 
time. 

As far as a free particle is concerned, the Langevin equation (1) may be 
written in terms of the linear momentum P = P(t) = mdX/dt as 

dP P 

— = - y _ + L(t). (2) 
dt m 

On the other hand, in the absence of inertial force, i.e., md 2 X/dt 2 = 0, the free 
Brownian particle is described by the stochastic differential equation 

dX 1 

Tt =-Ut). (3) 

In the Gaussian approximation it is assumed that the Langevin force L(t) 
present in (1 — 3) has the following statistical properties [2] 

(L(t)> = 0, (4) 

{L{t')L(t")) = D8{t -t"), (5) 

where D is a constant and the average value (... ) is evaluated from the probability 
distribution function associated with the environmental random function L(t). 
Equation (4) characterizes the irregularity feature of the Langevin force and Eq. 
(5) means that the interaction between the tagged particle and a generic 
environment is deemed to be Markovian in the sense that the autocorrelation 
function (5) is delta-correlated (white noise). 

For thermal open systems characterized by the temperature T and the 
Boltzmann constant k B , the Ornstein-Uhlenbeck process (2), along with (4) and 
(5), leads to the mean square momentum [3] 



AP(t) = 
N 



mk B T(l -e~Sr £ ), (6) 
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whereas the Fokker — Planck equation (the so-called diffusion equation) generated 
by the stochastic differential equation (3) with the statistical properties (4) and 
(5) yields the Einstein's mean square displacement [4] 



2k B T 

-?-t. (7) 



AX(t) = 

N 

Both Markovian results (6) and (7) are non-differentiable functions at t = 0. 
Hence, it has been claimed that physically there is no force applied to a free 
Brownian particle in the presence of inertial forces, as well as no velocity in the 
absence of inertial forces [5,6]. This fact suggests that the concept of trajectory of a 
free Brownian particle seems to be an elusive feature in the Markovian regime. 

The purpose of the present paper is to show how non-Markovian effects 
upon a free Brownian particle are responsible for the differentiability property of 
the mean square displacement (6) and momentum (7), thereby predicting an 
increase of the values of these physically observable quantities. 

Our paper discusses the Brownian motion of a free particle reckoning with 
two physical situations described by Fokker — Planck equations: In Sect. 2 the case 
of a Brownian particle in the presence of inertial forces (Ornstein — Uhlenbeck 
stochastic process) described by the non-Markovian Rayleigh equation is 
examined, whereas in Sect. 3 we solve our non-Markovian Smoluchowki equation 
describing a free Brownian particle in the absence of inertial effects. Concluding 
remarks are presented in Sect. 4. In addition, two appendices are included. 
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2 Brownian Motion in the Presence of Inertial Force 



A particle with mass m, position X = X(t) immersed in a generic environment 
(e.g., a non-thermal fluid) undergoes a jittering movement dubbed Brownian 
motion. The motion of this tagged particle may be mathematically described by the 
stochastic differential equation (1) (the so-called Langevin equation) in the form 

d 2 X dV(X) dX 

with the Langevin force given by L(t) = yftf^t), where the parameter b controls 
the environmental influence (fluctuations) upon the Brownian particle. Because 
the term bW{t) in (8) has dimensions of velocity, i.e., [length x time' 1 ], we can 
readily check that b may be expressed in dimensions of [length x time~ 1/2 ~\ and 
the function f(t) in dimensions of [time -1 / 2 ]. 

From the mathematical viewpoint the quantities X = X{t) and ¥ = f(t) in 
the Langevin equation (8) are interpreted as random variables belonging to the 
Kolmogorov probability space [7] in the sense that there exists a probability 
distribution function, T X w{x,\p,t), associated with the stochastic system {X, V], 
expressed in terms of the possible values x = {xj and \p = {i/^}, with i > 1, 
distributed about the sharp values q and cp of X and V, respectively. In addition, 
the average values of X and W are expressed as 

+oo 



(X)= jj xT xw (x,xp,t)dxdxp, (9) 



+ O0 



(W)= jj xpT xv (x,xp,t)dxdxp. (10) 

— 00 

For a free Brownian particle the Langevin equation (8) may be expressed in 
terms of the concept of linear momentum P = P(t) = mdX/dt as 



dP P 

-Tr = -Y- + Ybnt). (11) 
at m 



Noticing that 



t+s 



lim f (W(t'))dt' = 0, (12) 

£^0 J 



£ 

t 



with 



4 



+ 00 



<n0>= f xpT w (xp,t)dxp, 



the Fokker — Planck equation generated by (11) in the Gaussian approximation 
reads (see Appendix A) 



in terms of the probability distribution function g(p, t), the mean value of f(t) 

t+E 

<nO> = lim-[ <nO)dt', (14) 

£^0 £ J 

£ 

and the diffusion coefficient 

D p (t)=y£(t) (15) 

where the function £ (t), given by 

t+£ 



£(t) = —Kt) = — lim- <f(t )f(t ))dt dt , (16) 



t 



has dimensions of energy, i.e., [mass x length 2 x time -2 ]. Hence we call £ (t) the 
diffusion energy responsible for the Brownian motion of the particle dipped in a 
generic environment. 

An outstanding feature underlying the diffusion energy concept (16), which 
fulfils the validity condition 

< £(t) < oo, (17) 

is that it conveys in tandem fluctuation and dissipation phenomena through the 
function /(t) and the friction constant y, respectively. We can then state that (16) 
sets up a general fluctuation — dissipation relationship underlying all open systems 
described by the Langevin equation (11) and its corresponding Fokker — Planck 
equation (13). Both cases £(t) = and £(t) = oo are not concerned, for they may 
violate the validity condition of the fluctuation — dissipation relation (17). The 
former case may lead to dissipation without fluctuation, whilst the latter one may 
give rise to fluctuation without dissipation. 
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From the physical point of view the pivotal issue inherent in theory of 
Brownian motion is to determine the transport coefficient, that is, the time- 
dependent diffusion coefficient (15). At t = the answer to this question seems 
to be fairly straightforward, for the diffusion energy is null, i.e., £ (0) = y2)(0) = 0, 
meaning that there is no diffusive motion, D(0) = 0, associated with the initial 
condition g(p,t = 0) = S(p) to our Fokker — Planck equation (13). On the other 
hand, at the long-time regime, t -> oo, we assume that the steady diffusion energy 
can be identified with the thermal energy, k B T, at thermodynamic equilibrium, i.e., 

lim £{t) =^- = k B T, (18) 

t->co Z 

with the dimensionless function /(t) in (15) having the asymptotic behavior 

lim /(t) = 1. (19) 

t->co 

The physical significance of condition (19) has to do with the fact that 
environmental fluctuations decay to Markovian correlations, i.e., /(t) displays a 
local (short) range behavior in the steady regime. By contrast, non-Markovian 
effects show up at the range < t < oo. The correlational function 7(t)in (16) 
fulfilling condition (19) can be built up as (see Appendix B) 

/(t) = l-e"s (20) 

where the correlation time t c accounts for non-Markovian effects upon the 
Brownian particle. 

It is worth underscoring that the identity (18) expresses the principle of 
energy conservation, for the diffusion energy £ (oo) = yb 2 /2 comes from the 
Markovian Brownian dynamics (Langevin and Fokker — Planck equations at the 
steady regime), whereas the thermal energy k B T is a physical quantity stemming 
from the environment at thermodynamic equilibrium (equation of state for perfect 
gases: TV = Kk B T, where K is the number of particles within the volume V 
under the pressure T and temperature T; the constant k B denotes the Boltzmann 
constant and displays dimensions of energy per temperature). 

As a consequence of the conservation of energy (18), from (15) at the 
thermal equilibrium we can readily derive the Ornstein — Uhlenbeck diffusion 
constant as 

=Yk B T, (21) 
and the parameter b turns out to be of thermal nature, i.e., 

hk B T 

b = ± 

< Y 



(22) 
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Using (22), (20), (16), and (15) our Fokker — Planck equation (13) reads 



dgjy, t) 
dt 



^ Y dp ' 



— + 

m 



2k B T 



<no> 



) / -L\d 2 g(p,t) 

g(p,t)\ + Y k B T{l-e °^ - . (23) 



or 



^^ = m^ { ^ (p ' t)} + 7/CBr ( 1 - e 7^^' (24) 
in terms of the thermal position 



p = p + m 



2/c R r 



■<no>. 



(25) 



According to (23) a particle immersed in a heat bath undergoes a Brownian 
motion owing to the diffusion energy stemming from the thermal reservoir, k B T, 
multiplied by the correlational effects present in the function /(t) = (l — e~ t/£c ). 
In general, the non-equilibrium thermal energy £ (t) = k B T{l — e~ £,/£c ) is different 
from the equilibrium thermal energy, k B T. They equal only at the steady regime 
(18). In this thermal equilibrium state the stochastic process is called normal 
diffusion, whereas in the non-equilibrium thermal regime characterized by 
< /(t) < 1 the Brownian motion is said to be subdiffusive. 

It is worth pointing out that the diffusion constant (21) has been derived 
without specifying ab initio the form of the autocorrelation function of the 
Langevin force, 

(L(t')L(t")) = W(W(f)) = 2D p (oo)<nt')nt")>, 
as well as its average value, i.e., 

(L(t)> = ±^k B J{W{t)) = ± j2D p (oo)<<F(t)>. 

The crucial point is the long time behavior (19). This means that the arbitrariness 
as to the form of the statistical properties of the Langevin force L(t) leaves room to 
incorporate non-Markovian and averaging effects into the study of Brownian 
motion through the functions /(t) and (¥(t)) in the Fokker — Planck equation (23). 
So, starting from the initial condition g(p ' , t — 0) = S(p') a non-equilibrium 
solution to (23) reads 



,'2 



g(p', t) = 



(26) 



where 
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9(t) = 



mk B T 



1 

1 — e m — 



2yt c 



{2yt c - m) 



(27) 



The probability distribution function (26), which is expressed in terms of 
the evolution time t, the relaxation time t r = m/y, and the correlation time t c , 
gives rise to 



(P' 2 ) = 2Q(t)=mk B T 



<P') = o, 

1 — e m — 



2yt c 



(2yt c - m) 



e L c 



(28) 
(29) 



Making use of (29) the mean mechanical energy associated with the free 
Brownian particle is given by 



(P' 2 ) k B T 



(E(t)) = 



2m 



1 — e m — 



2yt c 



(2yt c - m) 



e f c 



(30) 



which reduces to 

k B T ( -iY t \ 
(E{t)) = -^-[l-e-^ 

in the Markovian limit. 

The mean square momentum, AP'(t) = ^(P' 2 ) — (P') 2 , reads 



(31) 



AP'(t) = AP(t) = 



mk B T 



1 — e m — 



(2yt c - m) 



(32) 



meaning that the average value of <F(t) in Eq. (25) has no influence upon the 
physically observable quantity (32). It is readily to check that the quantity 



dAP(t) 



dt 



2y 2 k B Tt c m - 2yt c 



t=o 



(2yt c — m) 2mk B Tyt c 



(33) 



does not diverge provided that < t c < m/2y, thereby implying that non- 
Markovian effects are responsible for enhancing the mean square momentum (32). 
Yet in the Markovian limit we find from (32) the result 



AP(t) = 



Jm/c B P ^1 — e m £ j 



(34) 



that is non-differentiable at t = 0. Non-Markovian effects in (32) therefore account 
for the existence of the concept of force acting upon a free Brownian particle in the 
presence of inertial force, since < t c < m/2y. 
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At the equilibrium regime, from (30) the energy equipartition is readily 
obtained as 

k B T 

(EM) = (35) 

while from the non-equilibrium solution (26) we derive the Maxwell — Boltzmann 
probability distribution function at thermal equilibrium 

g(p') = P ~2^7 ( 36 ) 

^j2nmk B T 

for the thermal momentum (25). 

3 Brownian Motion in the Absence of Inertial Force 

Now we start from the Langevin equation (8) in the Smoluchowski limit (the large 
friction case) 



P 

-Y — + YbV(t) 
m 





dP 


» 


dt 



(37) 



such that inertial effects in (8) can be negligible, i.e., md 2 X/dt 2 = 0. So the 
Brownian motion is approximated by the stochastic differential equation 

dX 1 dV(X) 

*=-rs i+6m (38) 

which gives rise to the following Fokker — Planck equation on configuration space 
in the Gaussian approximation (see Appendix A) 



df(x, t) _ d 
dt dx 



1 dV(x) 
y dx 



d z f(x t) 

f(x,t)+V x (t) \ (39) 



where the time-dependent diffusion coefficient is given by 

£{t) 

2>,(0 = -p (40) 

with the diffusion energy £ (t) being given by (16). Averaging effects of the random 
velocity, biW^t)), upon the drift coefficient bring about the time-dependent 
effective velocity 

1 dV(x) 

v eff (x,t) = -^+MnO>- (41) 

Y dx 

For thermal systems at thermodynamic equilibrium in which (18) is valid, 
the diffusion constant (40) becomes the Einstein relation 
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2), (oo) = 



k R T 



(42) 



Accordingly, the non-Markovian Fokker — Planck equation (39) in the presence of 
thermal fluctuations reads 



df (x, t) _ d 
dt dx 



1 dV{x) _ 
Y dx 



2k B T 



<no> 



f{x, t) + 1 {t) ■ 



(43) 



For (¥(t)) = an d /(t) = 1, the Fokker — Planck equation (43) yields the equation 
of motion 



df(x, t) _ 1 d 
dt y dx 



dV(x) 



dx 



f(x,t) 



k B Td 2 f{x,t) 
Y dx 2 



(44) 



early found out by Smoluchowski [8] and widely investigated in the literature [9 — 
14]. The Smoluchowski equation (44) is defined for Markovian correlations and no 
averaging effects in contrast to our Fokker — Planck equation (43) which takes into 
account non-Markovian effects in view of the time-dependent function l(t) as well 
as averaging effects present in the drift coefficient. 

Again, our account of the Einstein — Langevin stochastic approach presented 
above has pointed out that the main upshot of Einstein — the Founding Father of 
Brownian motion theory — rests on the fact that his diffusion coefficient (42) is a 
consequence of the conservation of energy (43) for thermal open systems at 
thermodynamic equilibrium. 

It is worth recalling that the Einstein relation (42) has been attained 
without specifying ab initio the form of the autocorrelation function of the 
Langevin force, 

k T 

<L(OL(o> = 2 Y 2 —(y<it')v<it")) = 2 Y 2 v x M(nt')nn), 



as well as its average value, i.e., 



k R T 



<L(t)> = ± Y 2-MnO> = ±X>/2ZUoo)<nt)>- 
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3.1 Non-Markovian Brownian Motion of a Free Particle 



As far as the non-inertial Brownian motion of a free particle is concerned, our non- 
Markovian Smoluchowski equation (43) with (20) at point x reads 



df(x',t) = _ 
dt + 



y dx y 



',t) k B T ( -L\d 2 f(x',t) 



Starting from the deterministic initial condition 

f(x',t = 0) = 5(x'), 



(46) 



characterized by £ (0) = and CF(0)) = 0, we obtain the following time-solution 
to (45) 



f{x,t) = 
in terms of the thermal position 

x = x + 



v yx 
e ^(0, 



AnA{t) 



2k R T 



and the function 



A(.t) = k.TjKW = k.T(t + tc .-t). 



Besides reducing to (46) at t = (along with t c = 0), solution (47) yields 

(X) = 0, 

and 



(X 2 ) = 



2k R T 



(t + t c e f c ), 



where the stochastic position X = X(i) corresponding to (48) is given by 



X = X (t) + 



2/c B r 



/ 



<<F(t))dt. 



Accordingly, the mean square displacement, AX(t) = 7 C^ 2 ) — W 2 < reads 



AX(t) = AX'(t) = 



2k B T 



(t + t c e £ c ). 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 
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We notice that averaging effects are unobservable, for they do have no influence 
upon the physically measurable quantity (52), albeit they can bring about a shift in 
the position (48). By contrast, non-Markovian effects account for enhancing the 
mean square displacement (52) of a free Brownian particle. 



For very short times t « t c , quantity (52) reduces to a constant value given 



by 



AX = 



2k B T 



(53) 



which is differentiable, i.e., dAX/dt = 0. On the other hand, in the Markovian limit, 
t c -> 0, quantity (52) reduces to the Einstein's renowned upshot 



AX(t) = 



2k B T 



(54) 



that is non-differentiable at t — 0, i.e., 



dAX(t) 



dt 



t=o 



k R T 



2yt 



— > 00. 



(55) 



t=o 



Hence, it is claimed that there is no concept of velocity of a Brownian particle in the 
strong friction regime [5,6]. In fact, the physically embarrassing outcome (55) is a 
consequence of neglecting non-Markovian features in (52). 
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4 Concluding Remarks 

In this paper we have examined the Brownian motion of a free particle immersed 
in a thermal reservoir reckoning with non-Markovian effects. In Sec. 2 we have 
shown that the mean square momentum (32) is a differentiable function in the 
presence of inertial force provided that non-Markovian effects increase the value 



In the absence of inertial force in Sec. 3 it has been predicted that non- 
Markovian effects enhance the mean square displacement (52), thereby assuring 
the mathematical property of differentiability of this physically observable 
quantity. 

It is worth point out that our chief results have been obtained without 
making use of the generalized Langevin equation 



put forward by Mori [15] on the basis of a close relationship between memory 
effects ingrained in the friction kernel /?(t' — t") and non-Markovian effects 
(colored noise) showing up in the autocorrelation function (L(t')L(t")), given by 



By contrast, our approach has predicted that non-Markovian effects 
independent of memory effects can be physically measured for a free Brownian 
particle immersed in a heat bath by means of the mean square displacement in the 
absence of inertial force as well as the mean square momentum in the presence of 
the inertial force. Surprisingly, this feature standing out in the Einstein — Langevin 
framework has been overlooked in the centenary literature about Brownian 
motion [1 — 20]. 



of (32). 




o 



{L(t')L{t"))=-(3(t -t"). 
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Appendix A: Generalized Fokker — Planck Equations 



Let us consider the stochastic differential equation 

dZ 1 dV(Z) 



dt 



Y dZ 



■ + aV(t). 



(Al) 



For Z = P, V(P) = Y 2 P 2 /2m, and a — yb Eq. (Al) reduces to the Langevin 
stochastic equation (11), whereas for Z = X, V = V(X), and a = b we obtain the 
Langevin equation (44). 



Equation (Al) gives rise to the Kolmogorov equation [7,11,17,18] 
df(z,t) 



dt 



= K/(z,t), 



(A2) 



where the Kolmogorovian operator K acts upon the probability distribution 
function f(z, t) according to 



Z(-l) fe 3" 
^-—^[A k {z,t)f{z,t)], 

k=l 



(A3) 



the coefficients A k (z, t) being given by 

[<(AZ) fe > 



A k {z,t) = lim 

£->0 



(A4) 



where the average values, ((AZ) fe ), are to be calculated about the sharp values z' 
from the probability distribution function 



T zw {z, xp,t)=S{z- z')T w {xp, t). 



(A5) 



According to the Pawula theorem [18], if the coefficients A k {z, t) in (A4) are 
finite for every k and if A k (z, t) = for some even k, then A k (z, t) = for all k > 3, 
thereby assuring the positivity of the probability density function f(z,t). So, if 
(z 2 - z 1 ) 3 ~0 such that (z 2 - z x ) 4 = 0, it follows then that A k {z, t) = 0, /c > 
3, where z 2 = z(t + £) and z x = z(t). Accordingly, the Kolmogorov equation (A2) 
can be approximated by the Fokker — Planck equation 



df(z,t) d Id 2 
J =~—M^> Of(z, t)] + -^[A 2 (z, t)f(z, t)], 



dt dz 
where the drift coefficient is given by 

[<AZ> 



A 1 (z, t) = lim 

£->0 



2dz< 



IdV 

= —— +a<no> 

Y dz 



(A6) 



(A7) 



and the diffusion coefficient reads 
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A 2 (z, t) = lim 

£->0 



<(AZ) 2 > 



t+e 



= a 2 \im^jj {V(t')V(t"))dt'dt' 

t 



with 



t+e 



lim - f 

£^o ej 



lim - <«F(t'))dr = <<f(0> 



(A8) 



(A9) 



and 



t+£ 

lim f {V(t'))dt' = 0. 

£^0 J 
t 



(A10) 



To evaluate the coefficients (A7) and (A8) we have used (Al) in the form 

t+£ 

edV f 

M = Z(t + £)-Z{t) = --— + a\ V(t')dt'. (All) 

Y dZ J 



For z= p, V(p) = Y 2 p 2 /2m, and a = yb the Fokker — Planck equation (A6) 
reduces to the Rayleigh equation (13), whereas for z= x, V = V(X), and a = b 
(A6) leads to the Smoluchowski equation (45). 
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Appendix B: The Time-dependent Diffusion Energy 



The time-dependent diffusion energy 

S(t) = £(oo)/(t) (Bl) 
in the Fokker — Planck equations (13) and (45) presents the correlational function 

t+E 

I(t) = lim-JJ (W(t')W(t"))dt'dt" . (B2) 

t 

On the condition that the autocorrelation function, (t')f(t")), can be 

given by 

/ (t'+t") \ 

(V(t')V(t")) = I 1 - e 2f c Jff(t'-t"), (B3) 

where t c is the correlation time of W(t) at times t' and t", it follows that (B2) 
becomes 

_t_ 

1(f) = 1 - e £ s (B4) 

reducing to /(oo) = 1 in the steady regime. 

The autocorrelation function (B3) defines a colored noise which in the 
Markovian limit, t c -> 0, changes into the so-called white noise 

(W(t')W(t")) = 8(t' - t"), (B6) 

meaning that the stochastic function W(t) is delta-correlated [16,19]. It has been 
argued that the Markov property (B6) is a highly idealized feature [19], because 
the physical interaction between the Brownian particle and the thermal bath 
always comes about for a finite correlation time. By the same token, a decade ago 
van Kampen has laconically stated: "Non-Markov is the rule, Markov is the 
exception" [20]. Hence, our auto-correlation function (B3) seems to be a more 
realistic feature of the Langevin force. 
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